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Abstract. Recently the concept of the pion cloud in the nucleon turned out 
to be successful in understanding the Gottfried sum rule violation observed 
by the New Muon Collaboration and the Drell-Yan asymmetry measured in 
NA51 at CERN. We propose a further possibility to test this concept at HERA 
through the analysis of the structure of deep inelastic scattering (DIS) events 
induced by pion-exchange. Momentum and energy distributions of outgoing 
nucleons as well as rapidity and multiplicity distributions are investigated us- 
ing Monte Carlo simulations. Most observables cannot distinguish this process 
from ordinary DIS, but in the energy distribution of final neutrons we find a 
significantly different prediction from the pion cloud model. Forward neutron 
calorimeters will be essential to test the concept of pions in the nucleon. 



1 Introduction 

In deep inelastic scattering (DIS) the incident lepton is scattered on a coloured quark. 
Normally this results in a colour field between the struck quark and the proton remnant, 
such that hadrons are produced in the whole rapidity region in between. In electron- 
proton collisions at HERA, this leads to particles being produced also close to the proton 
beam direction. The recent discovery by the ZEUS |TJ and HI || collaborations at HERA 
of large rapidity gap events has attracted much interest. This new class of DIS events 
have a large region of forward rapidity (i.e. close to the proton beam) where no particles 
or energy depositions are observed. The most forward hadronic activity being observed 
is then actually in the central part of the detectors. These large rapidity gap events 
cannot be described by standard models for DIS and hadronization || ||], Therefore 
the observation of a surprisingly large fraction (~ 10%) of events with a large rapidity 
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gap strongly suggests the presence of a final proton close to the beam momentum. These 
events have, therefore, been primarily interpreted in terms of pomeron exchange, although 
alternative models have recently been proposed 0. 

In this interpretation, the lepton interacts with a colorless object having the quantum 
numbers of the vacuum, i.e. the pomeron. The experimental signature is then a quasi- 
elastically scattered proton well separated in rapidity from the other produced particles. 
The leading proton escapes undetected by the main detector, but may be observed in 
leading proton spectrometers that are coming into operation in both ZEUS and HI. 

In the last few years, experiments on DIS have demonstrated that the internal struc- 
ture of the nucleon is more complicated than expected. The polarized DIS experiments 
performed by EMC and SMC at CERN have shown that only a small fraction of the 
proton spin is carried by the valence quarks ]7[. In addition, the strong violation of 
the Gottfried sum rule observed by NMC f8[ strongly suggests a d — u asymmetry of 
the nucleon sea. The new fits of the parton distributions M to the world deep inelastic 



and Drell-Yan data (including the dedicated NA51 experiment [ 10 1 ) seem to confirm the 
asymmetry. Both the violation of the Gottfried sum rule and the asymmetry measured 
in the Drell-Yan processes can be naturally accounted for by the presence of pions in the 



nucleon, as formulated in the pion cloud model [11 1. In view of these successes of the pion 



cloud model, it is mandatory to consider its role in other phenomena. 

The presence of such pions leads to an additional mechanism for nucleon production 
in DIS. In fixed target experiments, as (anti) neutrino deep inelastic scattering [12| for 



instance, it leads to the production of slow protons. The pion-exchange model describes 



the proton production on a neutron target [13| (extracted from deuteron data [Tj] obtained 
in bubble chamber experiments at CERN). With HERA kinematics the pion cloud induced 
mechanism leads to the production of rather fast forward protons and neutrons. The 
mechanism is shown schematically in Fig. 1. The virtual photon 'smashes' the virtual 
colorless pion into debris and the nucleon (proton or neutron) or an isobar is produced 
as a spectator of the reaction. In this respect there is full analogy to the reaction on the 
pomeron. Therefore, the pion cloud induced mechanism could also lead to rapidity gap 
events. In this paper, we investigate this processes and present quantitative results, not 
previously available in the literature. 

To understand these processes, not only protons but also neutrons in the forward di- 
rections are interesting |j!5| , [16| . Recently the ZEUS collaboration has installed a forward 
neutron calorimeter (FNCAL) |DJ which will provide additional experimental informa- 
tion. In analogy to the hadronic reaction pp — > nX, the pion-exchange is expected to be 
the dominant mechanism of the fast neutron production also at HERA |T5], [H]]. Thus, 
HERA open new possibilities to test the concept of pion exchange and the pionic structure 
in the nucleon. 

In the present paper we study several quantities which could be analyzed in HERA 
experiments; in particular using the main calorimeter, the leading proton spectrometer 



JT8[ (LPS) and the forward neutron calorimeter JT7| in ZEUS. The main aim of the study 
is to find the best signal to identify the discussed mechanism of scattering on a pion in 
the proton. 
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Figure 1: Fast forward nucleon production at HERA: (a) direct production through pion- 
exchange, (b) indirect production via a A resonance in pion- exchange, (c) pomeron- 
exchange. 

2 Pion— exchange mechanism of fast nucleon produc- 
tion 



In the meson cloud model jn| the nucleon is viewed as a quark core (called a bare nucleon) 
accompanied by the mesonic cloud. Restricting to the most important ttN component, 
the Fock state decomposition of the light-cone proton is 



|P> 



\(3q)) + / dyd 2 k T (j)(z,pT 



z,Pt) + 



mix 



z,Pt)) + 



with Z being the wave function renormalization constant which can be calculated by 
imposing the normalization condition (p\p) = 1. <J)(z,pt) is the light cone wave function 
of the ttN Fock state, where z is the longitudinal momentum fraction of the bare nucleon 
and Pt its transverse momentum. 

The presence of virtual pions in the nucleon leads to an additional mechanism for 
nucleon production referred to as 'direct spectator' (Fig. la) and 'sequential spectator' 
(Fig. lb) processes. The pion in the nucleon interacts with a virtual 7 producing a 
system X. For comparison we show the pomeron-exchange mechanism in Fig. lc. The 
cross section for the semi-inclusive spectator process ep — > e'NX can be written as 



dV p (ep -> e'NX) 



-/ttAt( 
Z 



z,t) 



da ev (x/(l - z)) 



(2) 



dxdQ 2 dzdp\ z j^^\ > > d(x/{l — z)) dQ 2 

The presence of the 7rA Fock component in the proton leads to the production of a 
spectator A which decays into a pion and a nucleon. The one-pion exchange contribution 
to the inclusive cross section can be obtained by integrating over unmeasured quantities 

da e P(x,Q 2 ) _ yi-* f t( -°' z) ^ f n _ ^da ew (x/(l - z),Q 2 ) 
Jo 



dx dQ 2 



dt f nN (l - z,t) 



d(x/(l -z))dQ 2 



(3) 



where a ew is the cross section for the inclusive deep inelastic scattering of the electron 
from the virtual pion. In practical calculations the on-mass-shell e7r cross section can be 
used. 
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The probability density to find a meson with light-cone momentum fraction x n = 
(1 — z) and four-momentum squared t (or alternatively transverse momentum p\ = — t(l — 
x n ) — m 2 x 2 ) is referred to as the splitting function, which quantifies the presence of virtual 
mesons in the nucleon. The splitting function f(x n ,t) to the tcN Fock state (Fig. la) is 



firN{ x iTy t) 



16tt 2 



.1 r 



-t)\F nN (x n ,t)\< 
(t — m 



2\2 



(4) 



and to the 7rA Fock state (Fig. lb) is 
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where M + = m/^+rriN and M_ 
via the isospin relations g 2 ^^ 
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■rriN ■ The couplings g 2 depend on the process, but 



g 



2 : 1 and ^ 7r+A o : g*o A+ : g\ 



-A+^ 



2 : 3 there are only two independent couplings which we take as 9n—MrOt> 

/4?r = 13.6 pi] and 



g 



-A++ 



/4tt = 12.3 GeV- 2 



The -PmbOetd are vertex form factors, which account 
for the extended nature of the hadrons involved. The form factors used in meson exchange 
models are usually taken to be functions of t only. As discussed in ref. [[B| such form 
factors are a source of momentum sum rule violation and it was therefore suggested to use 
form factors which are functions of the invariant mass of the intermediate meson-baryon 



system, i.e. M 2 MB (x^p 2 T ) = ^ + 



It can be shown that such a vertex function arises naturally if one computes the 
splitting function f(x n ,t) in time-ordered perturbation theory in the infinite momentum 
frame |22| . This functional form is typical for parameterizing the light-cone wave function 
of composed systems (see e.g. p3|). 



In all calculations discussed below the vertex form factors have been assumed in the 
exponential form 



F M B{x-n,p 2 T ) = exp 



9A 2 



By using the kinematical relation [24 
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(7) 



the form factor given by Eq.(|G]) can be equivalently expressed in terms of x n and t in the 
simple form: 



F M B(x n ,t) = exp 



m; 



The cut-off parameters used in the present calculation (A nN = 1.10 GeV and A nA = 
0.98 GeV) have been determined from the analysis of the particle spectra for high-energy 
neutron and A production [jID|, i.e. pp — > nX and pp — > A ++ X. With these cut-off 
parameters the NMC result for the Gottfried sum rule || which depends sensitively on 
Amb, has been reproduced [T9[]. Furthermore the model describes the u — d asymmetry 
extracted recently from the Drell-Yan experiment NA51 at CERN ||25|| . We note, however, 
that all results of this paper would be quite similar if traditional dipole form factors with 
cut-off parameter of 1.0-1.2 GeV had been used instead of Eq. (fl). 
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In hadronic reactions quite often the Regge approach was used rather than the light- 
cone approach. In order to obtain the flux factor in the Regge approach it is sufficient to 
replace in Eq.(^) x T by x\~ 2oilT( - t \ where the pion's Regge trajectory a n (t) = a n (t — m 2 ). 
The reggeization is important for small i T and/or large t. This is a kinematical region 
where the flux factor, especially with the vertex form factor Eq. (||), is rather small. 
Furthermore in the Regge approach, in contrast to the light-cone approach, it is not clear 
whether it would be fully consistent in the lepton DIS to use the on-shell pion structure 
function. However, since the difference is important only in very limited region of the 
phase space, in practice both approches lead to almost identical flux factors. 



3 Results and Discussion 

The formalism presented above has been implemented in the Monte Carlo program 
Pompyt version 2.3 f26f . This program, which was originally for diffractive interac- 
tions via pomeron exchange, simulates the interaction dynamics resulting in the complete 
final state of particles. The basic hard scattering and perturbative QCD parton emission 
processes are treated based on the program Pythia |27) and the subsequent hadroniza- 
tion is according to the Lund string model @] in its Monte Carlo implementation Jetset 



f2"?fl which also handles particle decays. 



The main difference in comparison to the pomeron case is the replacement of the 
pomeron flux factor by the pion flux factors given by Eqs. (§,0) and the pomeron structure 
function by the pion structure function. The pion case is better contrained than the 
pomeron case, due to the better known pion structure function where those for the on- 
shell pion can be used. The pion parton densities from the parametrisation GRV-P HO 
(MS) is therefore used. It is important to mention in this context that the absolute 
normalization of the cross section for the production of the spectator nucleon via pion- 
exchange mechanism depends on the absolute value of the pion structure function. At 
the small- x relevant at HERA, the structure function is completely dominated by the 
pion sea contribution which is not very well known. Experimentally the pion structure 



function can be determined from the Drell-Yan processes only for x > 0.1 [^, pOfl . If 
the pion-exchange mechanism is the dominant mechanism of fast neutron production, 
the coincidence measurement of scattered electrons and forward neutrons may allow the 



determination of the pion deep inelastic structure function ||16|| . When considering the 
event structure, however, the precise value of F% is not required. 

When the deep inelastic scattering is on a valence quark (antiquark) the pion rem- 
nant is simply the remaining antiquark (quark). A colour triplet string is then stretched 
between them and hadronization described with the Lund string model 0]. In case it 
is a sea quark (antiquark) that was struck, the pion remnant contains the associated 
sea antiquark (quark) in addition to the valence quark and antiquark. A string is then 
stretched between the struck quark (antiquark) and a valence antiquark (quark), whereas 
the remaining valence quark (antiquark) forms a meson together with the spectator sea 
antiquark (quark). 

For the results presented below we have made simulations corresponding to the HERA 
conditions, i.e. 26.7 GeV electrons on 820 GeV protons. The results for the above pion 
exchange mechanism are compared with normal DIS on the proton, which is simulated 



5 



protons 
ontiprotons 




neutrons 
ontineutrons 




b) 



spectator neutrons 
(7T + - exchange) 
neutrons from 
DIS on n + 
antineutrons from 
DIS on 7T + 



c) 



200 400 600 8000 
E p (GeV) 



200 400 600 8000 
E n (GeV) 



200 400 600 800 
E n (GeV) 



Figure 2: Energy spectra in the HERA lab frame for nucleons (p, p,n,n) from (a,b) conven- 
tional DIS on a proton (obtained with Lepto^ and (c) DIS on an exchanged tt + (obtained 
with POMPYT ). 

with Lepto 6.3 fH using the MRS(D-') parton distributions W^ - in all cases, events 



are simulated according to the cross section formulae and are constrained to be in the 
kinematical region x > 10~ 5 , Q 2 > 4 GeV 2 . 

In Fig. 2 we show the resulting energy spectra of nucleons (p, p, n, n) in the lab frame 
of HERA. This is of direct interest for measurements in the leading proton spectrometer 



18| and forward neutron calorimeter |17|]. Neutrons from the pion exchange mechanism 
have large energies giving a spectrum with a broad peak around E « 0.7E beam , i.e. 
around 500 GeV, whereas the corresponding spectrum from DIS on the proton decreases 
monotonically with increasing neutron energy. In the region of interest, say 400-700 GeV, 
the two processes have a similar absolute magnitude. An observable effect from DIS on a 
pion should be therefore possible. 

While the energy distribution of primary A's is very similar to that of the direct neu- 
tron production after the A — > nn decay the energy distribution of the secondary 



nucleons becomes peaked at smaller energies of about 400 GeV ||32|| . The two-step mech- 
anism is, however, much less important for the production of neutrons. First of all the 
probability of the 7rA Fock states in the light-cone nucleon wave function is much smaller 
than the probability of the nN component: ~ 0.3P n N jl^D - Secondly, the isospin 

Clebsch-Gordan coefficients favour the decay of the A into the proton over the decay 
into the neutron channel with the proton/neutron branching ratio | : |. The analogous 
branching ratio for the direct component is | : |. All this imply that both the 1-step and 
2-step mechanisms produce comparable amounts of protons. In contrast, the two-step 
mechanism produces about 10 times less neutrons than the 1-step mechanism. This means 
that in a first approximation the two-step process may be neglected for the spectrum of 
neutrons. Therefore we concentrate on the comparison of DIS on tt + , having a neutron 
as spectator, with standard DIS on the proton. 

The calculated transverse momentum (p?) distributions are shown in Fig. 3. As can 
be seen, the distribution of the spectator neutrons falls faster with increasing p\ than 
that from standard DIS. It can be expected that the distribution of neutrons produced in 
the two-step process in Fig. lb is less steep than those produced in the direct process in 



6 



> 

O 
\ 

c 



Q_ 
X> 
\ 

b 



10' 



10- 



10 



spectator neutrons 
(n* - exchange) 

neutrons from hadronization 
(LEPTO 6.3) 




0.2 0.4 0.6 0.8 



1.2 1.4 1.6 



1.8 2 
Pt 2 (GeV 2 ) 



Figure 3: Distribution in transverse momentum of neutrons from DIS on the proton 
(dashed histogram from Lepto^) and from DIS on a n + with a neutron as a spectator 
(solid histogram from Pompyt^). 



Fig. la. The higher overall level of DIS on the proton can be reduced by a cut in neutron 
energy, as is obvious from Fig. 2. Still, the difference in shape of the p^-spectra in Fig. 3 
is presumably too small to be exploited experimentally. A safe conclusion does, however, 
require further analysis including, e.g., finite angular acceptance of FNCAL. 

To study other characteristics of events arising through DIS on a virtual pion and 
compare with standard DIS on the proton, we consider spectra of different quantities 
normalized as 

rt , 1 dN 

where k can be any kinematical variable and N event is the number of events. This give 
emphasis to shapes irrespectively of normalisation and statistics (of data and Monte Carlo 
samples) . 

A quantity with especially nice transformation properties under longitudinal boosts is 
rapidity defined as 

' /':. \ 



y 




(10) 
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Figure 4: Rapidity distributions of (a) all stable particles, charged pions and j's and (b) 
protons and neutrons from DIS on the proton as obtained from Lepto. 

where E is the energy and p z the longitudinal momentum along the proton beam axis. 
For massless particles this quantity is identical to the pseudo-rapidity defined by 

rj = -lntan(9 /2) (11) 

where 9 is the angle of a particle with respect to the proton beam, i.e. rj > is the proton 
hemisphere in the HERA lab frame. 

In Fig. 4 and 5 we show the pseudo-rapidity distributions of different particle species 
produced in DIS on the proton and DIS on a 7r + , respectively. In Fig. 5a spectator 
neutrons are not included, but shown separately in Fig. 5b. For example, the size of 
the beam pipe hole in FCAL (9 = 1.5°), assures that in almost 100% of the specta- 
tor nucleons (proton/neutron) leaves the main ZEUS detector without any energy loss. 
As seen by comparing Fig. 4a and Fig. 5a the pseudo-rapidity spectra of n ± and 7 are 
rather similar in the two cases. The pseudo-rapidity spectrum of spectator neutrons 
(Fig. 5b) has a maximum at only a slightly higher value compared to the peak of neutrons 
from non-diffractive DIS on the proton (Fig. 4b). These predicted neutron distributions 
should be considered in the context of the pseudo-rapidity coverage of the forward neu- 
tron calorimeter. In general, the neutron acceptance is a complicated function of both 
polar and azimuthal angle. The ZEUS FNCAL geometry limits pseudo-rapidity coverage 
approximately to 7 < r] < 10. The Lund hadronization model predicts a small amount of 
nucleon-antinucleon pairs produced in DIS on the pion (Fig.5cd). 

The pseudo-rapidity variable is of particular interest in the context of large rapidity 
gap events. These have been defined by i] max giving, in each event, the maximum pseudo- 
rapidity where an energy deposition is observed. Based on our Monte Carlo simulated 
events using Lepto and Pompyt we extract this f] ma x- variable and show its distribution 
in Fig. 6 for conventional non-diffractive DIS on the proton, DIS on an exchanged 7r + and 
diffractive DIS on a pomeron. Since our aim here is to demonstrate the genuine physics 
effects of the models, we have not included any experimental acceptance effects or rapidity 
gap requirements in this study. Doing this will severely distort the distributions at large 
r] max and, therefore, one cannot make direct comparisons with the available measured 
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Figure 5: Rapidity distributions of the specified particles produced in DIS on a ir + (neutron 
spectator) as obtained from Pompyt. 

distributions. Thus, from this model study we find a shift of about one unit towards 
smaller f] max i n case of DIS on the pion as compared to normal DIS. For f] max < 6, these 
two processes contribute about equally to the rate. 

For the spectrum of rj max for DIS on the pomeron, we have taken a set of parameters 
which is usually called 'hard pomeron' in the literature. The pomeron is assumed to 
contain equal amounts of the light quarks, i.e. u = u = d = d, each with a density 
distribution 

zq(z) = -z(l - z) , (12) 

with the normalization chosen such that the parton distributions fulfill the momentum 
sum rule. The pomeron flux factor is here taken as the ratio of the single diffractive cross 
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Figure 6: Distribution of r] max (see text) in non-diffractive DIS on the proton (solid), in 
DIS on the virtual 7r + (dashed) and in DIS on the pomeron (dotted); pure physics of the 
models without experimental acceptance effects. 



section and the pomeron-proton total cross section [ 33 1 



where the simple parameterization is obtained by fitting the numerator to single diffractive 
cross section and the denominator is taken as a(JPp — > X) = 2.3 mb obtained from a Regge 
analysis p4j of elastic and single diffractive scattering. The resulting 7] max distribution 



from DIS on the pomeron is considerably different from the other two cases. 

From Fig. 6 one may conclude that the pion-exchange induced DIS leads to events 
with intermediate size rapidity gaps rather than to those with large gaps. Nonetheless it 
is important to verify experimentally the effect of the pion exchange by, e.g., correlating 
Vmax with fast forward neutrons measured in FNCAL (for technical details see [|35|j). 

The flux factor given by Eq. (H) with a cut-off parameter of the vertex form factor 



extracted from the high-energy neutron production data |L9j predicts that the pion carries, 
on average, a fraction 0.3 of the proton beam momentum |[22|| . This implies that as a first 
approximation the pion-induced DIS processes can be viewed as an electron scattering on 
the pion with effective energy E e ff ps 0.3 • E beam . Because of the smaller energy of the 
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Figure 7: Multiplicity distributions of (a) all stable particles, (b) charged pions from DIS 
on a 7r + (full curves) and DIS on a proton (dashed curves). In (c) j's from DIS on the 
proton and (d) 7's from DIS on n + . 



pion one could expect smaller multiplicity of electron-pion DIS events in comparison to 
those for the electron-proton DIS. In Fig. 7 we compare the model predictions (without 
experimental acceptance effects) of the multiplicity spectra for DIS on the proton with 
those on ir + . 

The multiplicities in DIS events on the pion is noticeably smaller than in DIS events 
on the proton; the average multiplicity is about 20 and 30, respectively. The dominant 
contribution to the multiplicity spectra comes from charged pions (11.6 on the proton vs. 
7.5 on the 7r + ) and 7's (14.1 on the proton vs. 8.0 on 7r + ). The even-odd fluctuations 
of the multiplicity spectra of photons is not statistical, but caused mainly by the decay 
7t° — > 77. Thus, as expected the multiplicity of standard DIS events is typically larger 
than in pion-induced DIS events. However, due to the large fluctuations in multiplicity 
and the overlap between the distributions for the two cases, as well as the distortions that 
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limited experimental acceptance will create, it is not clear whether this difference can be 
used as a discriminator. This needs further considerations. 



4 Conclusions 



in 



The concept of a pion cloud in the nucleon was recently found to be very useful [JT9J, |25 
understanding the Gottfried sum rule violation observed by the New Muon Collaboration 
H and the Drell-Yan asymmetry measured recently in the NA51 Drell-Yan experiment 



at CERN [llj. In the present paper we have investigated several quantities in order 



to find useful observables which would help to verify this concept using deep inelastic 
electron-proton scattering at HERA. We have therefore analyzed the structure of deep 
inelastic events induced by the pion-exchange mechanism. In particular, we have studied 
distributions of final nucleons as well as rapidity and multiplicity spectra. 

Most of the event characteristics do not provide a direct possibility to distinguish 
the events from DIS on a pion from the ordinary events with DIS on a proton. A clear 
difference is, however, found in the energy spectrum of outgoing neutrons. We find that the 
pion cloud model predicts an energy distribution of neutrons which substantially differs 
from the standard hadronization models. While the pion-exchange mechanism leads to an 
energy spectrum which peaks at an energy of about 0.7E beam , i.e. at about 500 GeV, the 
spectrum of neutrons produced in the standard hadronization process following DIS on 
the proton decreases monotonically with increasing neutron energy. This should facilitate 
to discriminate between the two processes, in particular since they have cross sections of 
similar magnitude in this energy region. Therefore, the experiments with forward neutron 
calorimeters should shed new light on the nucleon structure in terms of a pion content. 

We have shown that the pion cloud induced mechanism practically does not contribute 
to the large rapidity gap events observed recently by the ZEUS and HI collaborations [I], |2j 
and cannot be a severly competing mechanism for the pomeron exchange. The multiplicity 
of the pion cloud induced events is about 60-70% of that for standard hadronization on 
the proton, but given the large fluctuations it is not clear to what extent this difference 
can be exploited. 

Our results on the pion exchange mechanism are more general than the detailed for- 
mulation of the pion cloud model. Since essentially the same pion flux factor is obtained 
in Regge phenomenology, our results may also be taken as a representation of Regge-based 
expectations. 

In this study we have omitted experimental effects due to finite appertures, clustering 
effects in the main detector, finite energy thresholds, detector efficiencies, etc., which may 
distort the observed spectra. Many of them are quite important in order to understand 



and interpret the observed spectra and we plan a future study [36] of such effects. 
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